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THE EFFECT OF TORSIONAL FLEXIBILITY ON THE ROLLING 
CHARACTERISTICS AT SUPERSONIC SPEEDS 
OF TAPERED UNSWEPT WINGS 
By Warren A. Tucker and Robert L- Nelson 

SUMMARY 


An analysis is presented of the effect of torsional flexibility on 
the rolling characteristics at supersonic speeds of tapered unswept 
wings with partial-span constant -percent -chord ailerons extending inboard 
from the wing tip. The geometric variables considered are aspect ratio, 
taper ratio, aileron span, and aileron chord. 

The shape of the wing-torsional -stiffness curve is assumed and the 
twisting moment is considered to result solely from the pressure distri- 
bution caused by aileron deflection, so that the necessity of using a 
successive -approximation method is avoided. 

Because of the complexity of the equations resulting from the 
analysis, numerical calculations from the equations are presented in a 
series of figures. A computational form is provided to be used in con- 
junction with these figures so that calculations can be made without 
reference to the analysis. 


iNTRODUCTION 


In a previous paper (reference l) an analysis was presented of the 
effect of torsional flexibility on the rolling characteristics at super- 
sonic speeds of rectangular wing-aileron combinations. The present 
paper deals with the more general case of linearly tapered wings having 
straight midchord lines and incorporating partial-span constant -percent - 
chord ailerons extending inboard from the wing tip (see fig. l) and may 
thus be considered a companion piece to reference ,2, which gives an 
analysis of the rolling characteristics of this class of wings at sub- 
sonic speeds. 

The necessity of using successive approximations is avoided by 
assuming, as in reference 1, that the principal twisting moment is 
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contributed by the loads arising from aileron deflection. The various 
aerodynamic forces and moments are calculated by methods based mainly 
on the linearized equations for three -dimensional supersonic flow- 

The geometric variables considered are aspect ratio, taper ratio, 
aileron span, and aileron chord. The torsional stiffness of the wing 
is assumed to vary inversely with the cube of the distance outboard from 
the wing center line- 

Because of the length of the equations for the various quantities 
appearing in the analysis, numerical calculations have been made for 
ranges of the variables considered. The results are presented in a series 
of figures and can be used, together with a computational form that is 
provided, to make calculations without reference to the details of the 
analysis. 
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SYMBOLS 


aspect ratio 
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wing span 

total span of two ailerons 

wing chord at spanwise station y 

aileron chord at spanwise station y 

root chord 

tip chord 

taper ratio (c^j 


slope of hinge line (see fig. 1) 



slope of leading edge (see fig. 1) 



wing area 
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Bemlspan (b/2) 
coordinate axes (see fig* l) 


free -stream Mach number 

free-stream Mach number component normal to 
hinge line 


ratio of specific heat at constant pressure to 
specific heat at constant volume (1.40 for air) 

free-stream static pressure 

free-stream dynamic pressure ^|-p 0 M^ 
lifting pressure 

lifting-pressure coefficient (p/q) 
additional two-dimensional lifting-pressure 


coefficient due to aileron deflection 


48 


3n 


1 - r 

deflection of one aileron, measured in plane 
perpendicular to wing and parallel to free 
stream, positive when right aileron is 
deflected up 

angle of twist at any spanwise station, positive 
when leading edge of right wing moves up 

angle of twist at reference station (midspan of 
aileron) 
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T' = 
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GJ 


wing-tip helix angle 

rolling angular velocity, positive when right 
wing moves down 

constant of proportionality for parabolic twist 
distribution 

vertical disturbance velocity 
free -stream velocity 

rolling moment, positive when tending to depress 
right wing 

rolling-moment coefficient (L/qSb) 
rolling-mcment-effectiveness coefficient for two 
ailerons 



damping-in-roll coefficient 



/°C 

rolling-moment -loss coefficient r-r* 
. • \ d0 r 


aerodynamic twisting moment at any spanwise 
is 

m dyl 


station 



pitching moment about midchord at any spanwise 
station 

torsional-stiffness parameter 



“0 


torsional-stiffness parameter (t/©) 
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Subscripts: 

E 

0 

1 , 2 , 3 , b 

"bj c 


torsional-stiffness parameter at reference station 
(midspan of aileron) 

constant of proportionality for stiffness distri- 
bution 

trailing-edge-angle reduction factor 

included trailing -edge angle measured parallel 
to free Stream 

included trailing -edge angle measured normal to 
hinge line 


rigid 
sea level 

used to indicate various regions of the wing 
used to indicate various stations 


All angles are in radians, unless otherwise specified. 

ANALYSIS 

General 


In a steady roll, the resultant rolling moment acting on a wing of 
any plan form is zero. This statement can be written in equation form as 


+ PC ‘p | + fit V r - 0 


( 1 ) 


Equation (l) can be solved either for the rate of roll per unit 
aileron deflection or for the ratio of the rate of roll of a flexible 
wing to the rate of roll of a perfectly rigid wing (d r = 0) for a 
given aileron deflection. The following equations are obtained: 
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The problem now, becomes one of finding expressions for the quantities 0C^ 
PCj , 0Cj , and 0 r /S for the class of wing-aileron combination 

.r 

considered; the rest of the analysis is devoted to this task. 

A knowledge of the twist along the span is needed for the determina- 
tion of PCi and 0 r /8; this requirement in turn necessitates finding 

e T 

the aerodynamic twisting moment at any spanwise station. In order to 
evaluate the aerodynamic twisting moment the following assumptions are 
made: 

(1) The elastic axis and the center of pressure at any spanwise 
station of the loads resulting frcrn rolling and from twist lie very near 
the midchord line. This assumption allows the twisting moments resulting 
from rolling and from twist to be neglected. 

(2) The aerodynamic twisting moment results solely from the loads 
caused by aileron deflection. This assumption is a direct result of the 
first one. 

( 3 ) The loads caused by aileron deflection, wing twist, and rolling 
are mutually independent.. This assumption corresponds to the familiar 
superposition principle of linearized theory. 

(4) The aileron does not twist with respect to the wing. 

An important result of making these assumptions is the elimination 
of any necessity for using a method of successive approximations. 

The main part of the analysis follows: First, the pertinent pressure 

distributions are discussed; next, the method of finding the twist along 
the span is outlined; and last, the results for the various quantities 
are given. Since the mathematics involved in each step is lengthy, no 
attempt is made to give all the details. Those readers whose primary 
interest is in the application of the results are advised to turn to the 
section entitled "Computational Procedure and Discussion" and to table I. 


Pressure Distributions 


Three pressure distributions are necessary in order to determine 
0%, 0 c lp> ^<1 & c Iq * These pressure distributions are, respectively, 
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that due to aileron deflection, that due to a wing twist increasing 
linearly with distance along the span (corresponding to a rolling wing), 
and that due to a wing twist increasing with the square of the distance 
along the span (the reason for using a parabolic twist is discussed in 
the section entitled "Twist ratio"). Each pressure distribution can be 
found by use of the methods of references 3 and h if appropriate source - 
sink distributions are used to represent in one case the deflected / 

aileron and in the other cases the twisted wing. The details of the 
derivations are not given} an illustration of the use of the method for 
finding the pressure distribution over a twisted wing can be found in the 
appendix of reference 1. 

The equations for the pressure distribution due to aileron 
deflection are “3 



for the regions indicated in the following sketch: 



Sketch 1 
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The equations for the pressure distribution due to a twist increasing 
linearly with distance along the span (w = py) are 



x - k 2 s 


and the regions are indicated in the following sketch: 



x 


Sketch 2 
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The equations for the pressure distribution due to a twist increasing 
with the square of the distance along the span (w = ay^) are 



x - kgS 

and the regions are the same as for the wing with twist increasing linearly 
with distance along the span (see sketch 2). 
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Twist 

Aerodynamic twisting moment . - The aerodynamic twisting moment T at 
any spanwise station y is simply the total moment acting outboard of 
the station. Thus, 



where m, the pitching moment about the midchord line at any spanwise 
station, arises from the pressure distribution caused by aileron 
deflection and can be calculated from equations (h) • 

The resulting expressions for the aerodynamic twisting moment were 
extremely complicated. The twisting -moment distribution was therefore 
approximated by a series of straight lines as shown in figure 2, where 
the exact and approximate distributions for a particular configuration 
are compared for two values of 0. The equations for T' a , T ' , and T' 

(see fig. 2) are as follows: 
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The spanwise stations for T' a , T'^, and T' c are given by the 
following equations: 




_2X Ca _ 1 - g ha 

1 - \ 1 - | c . 1 - | h 




(9) 


The straight-line approximations to the twisting-moment distribution 
(see fig. 2) are given by 
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Twist at any spanwise station .- The twist 9 at any spanwise 
station y can be calculated from the equation 


9 


py 
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d0 
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dy 


(n) 


where dfl/dy can be found from the definition of the torsional-stiffness 
parameter GJ, which is 


T 

dfl/dy 


GJ = 


( 12 ) 
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and is related to the torsional -stiffness parameter mg by the following 
equation: 


1_ 

“0 


py 


u 


GJ 


(13) 


If the assumption is made that 


± p 

“ = ky 2 


GJ 


(14a) 


then 


i = K3 
“e J 


(l4b) 


This variation of torsional stiffness is known to be typical of many 
subsonic fighter airplanes (reference 2) and may be expected to approxi- 
mate the stiffness distribution of airplanes or missiles designed for 
supersonic flight. 

If equation (l4a) is substituted in equation (12) and the result 
substituted in equation (ll ) t then the twist is expressed as 


0(j) = k / y 2 T(y) dy 


(15) 


If a change is made in the variable of integration by writing 


y 

=7 


and if T 1 is defined as 


T ' = T 


23 




( 16 ) 


then the expression for 9 becomes 


C kqs3 c . 3 P n 


eOi) = 


23 
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(17) 
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The integration must of course he carried out in steps, since T' is not 
defined continuously over the span. 

The subscripts 1, 2, 3, and 4 may conveniently be used to 
indicate the regions designated by the same subscripts in figure 2. 

After substitution of the values of T 1 from equations (10) in equation (17) 
the twist equations become: 
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(18a) 

(18b) 



9 = T 'a A 3 _ tA / t * b - T * a \ 

^ C p Jcqs 3 c t 3 1 " H a \3 4/ \ \ " \ ) 12 


(18c) 



(l8d) 


NACA TN 1890 


15 


Twist ratio .- For purposes of calculation, the twist may he expressed 
as the ratio of the twist 9 at any station to the twist 9 r at a 


reference station, taken at the mldspari of the aileron 


= l-±^ 

2 b 


The twist is plotted in this manner in figure 3 for a particular configu- 
ration for two values of 3* Because 0/0 T is not defined continuously 
over the span and because the defining expression is very complicated in 
each region, the subsequent evaluation of 0C is greatly simplified 

if a simple approximate equation for 6 /Q T) continuous over the span, is 
used. This approximation, which is also plotted in figure 3, is given 
by the following equation: 



. 1 

2 


1 - 


V 2 


(19) 


Derivatives 

EoUing-moment-effectiveness coefficient 0C i& .- The rolling-moment- 
effectiveness coefficient pCjg is first found for a wing-aileron con- 
figuration of zero thickness (or zero trailing- edge angle) by suitable 
integration of the pressure distribution given by equations (4) . A 
factor F is then applied to account for the effect of the trailing- 
edge angle . 



result obtained for the wing of zero thickness is given by the following equation: 
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Damping- in-roll coefficient PC^.- The damping- in- roll coeffi- 
cient PC 2 p is obtained "by proper integration of the pressure distribu- 
tion given by equations ( 5 ) • The following result is obtained: 


-pc 


Zp 


2 

jt(l + X)(l - X) 3 


l-i-itiLt ^. 4 003-15 + 35 - i °£ 3 - 8 ^ 
3(1 - I 2 ) 7 ' 2 9(1 - I 2 ) 3 


X 1 * 2 - 64 + 55 2 _ 1-2$ 

3(1 + X)(l - X)3 2(1 - O 3 3(1 + X)(l - X) 2 (1 - |) 2 


2X 2 1 

1 - X 2 1 - I 


( 21 ) 


Roll tng-moment- loss coefficient pC iq .- The rolling-mcment- loss 


is obtained by suitable integration of the pressure 


coefficient PC 2 ( 
distribution giveji by equations (6) by using 



( 22 ) 


from equation (l 9 ), since 9 - The following expression results: 



1 ^a'j 2 2 

2 V 15 * (l - X) 4 (l + X) 


3(1 - t 2 ) 4 <9 ’ 361 
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j2. 6 - 24£ + 3‘3l 2 - 20^3 

30(1 - X) 4 (l + X) (1 . jit 

X 4 6 - 18 1 + 17..I 2 . V. 2X 3 1 - 2j 

6(1 - x) 3 (i + x) (1 - I) 3 (1 - X) 2 (1 + X) (1 - if 


■ 2X 2 1_ 

(1 - X)(l + X) 1 - I 


(23) 
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Twist parameter 9 r /b.~ The twist parameter 0 r /5 can he evaluated 


from equations (l8c) and (l8d) for 1 - i less than or greater 

than T) a , respectively. In order to perform the evaluation, the value 
of Cp is taken from equations (4) to he the following: 


45 


f 7 ? 


P 


(24) 


The use of this linearized-theory result is conservative. The value 
of k is found from equation (l4h) to he the following: 


k = 


24 


^ ‘ 2 ¥)' 


(25) 


The following expressions are obtained for 0 r /5: For 1 - A — < n 

A 1 2 *b a 9 
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(26a) 


and for 1 - | ^ > T) q , 
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' COMPUTATIONAL PROCEDURE AND DISCUSSION 


In order to simplify the solution of equations (2) and ( 3 ) a compu- 
tational form has "been prepared. This computational form, which is presented 
as table I, is practically self-explanatory and can "be used without refer- 
ence to the analysis • Figures ^ to 15 contain the necessary information 
to he entered in the table. Most of these figures present the results of 
numerical calculations made from equations (20), (21), ( 23 ), and (2.6) j 
the variables I and £ appearing in these equations have been replaced 
by the appropriate combinations of (3A, X, and c a /c so that the figures 
contain these more familiar quantities. The curves for the rectangular 
wing (X = 1.0) were taken from reference 1. Although in principle these 
curves could have been obtained from the equations of the present paper, 
the work involved in evaluating the many indeterminate forms arising 
when X = 1.0 made such a procedure impractical. 

The factor F which appears in equation (20) is introduced to 
allow for the effect of the trailing-edge angle. This factor is the same 
one shown, in figure 7 of reference 1 and was derived by using Busemann's 
second-order approximation for isentropic flow. For the purposes of the 
present paper, the Mach number and trailing-edge angle must be taken 
normal to the aileron hinge line. The factor F is plotted in figure 10 
against M^, the Mach number component normal to the aileron hinge line, 
with the included trailing-edge angle measured normal to the aileron 

hinge line, as the parameter. The relations between 1% and M and 
between $N and $ are given by the following equations: 



For ailerons other than those having flat sides, the trailing-edge angle 
should be taken as the angle between straight lines drawn from the 
trailing edge to the airfoil surface at the hinge line. For this 
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procedure can be shown to give the correct value of F independent of 
the airfoil profile. However, this method of determining F for curved- 
profile control surfaces is not valid for trailing- edge -angle corrections 
to hinge moments . 

The aileron deflection 6 'used in is measured in a plane 

parallel to the free stream. The value of the rolling-moment- effectiveness 
coefficient referred to an aileron deflection measured perpendicular to 
the hinge line can be obtained, if desired, by dividing 0C2g by the 
quantity \fl + k] 2 . 

The minimum value of 0A for which each curve is plotted corre- 
sponds in most cases to the limit of applicability of the equation for 
the curve. Extrapolation of the curves is not recommended. The maximum 
value of b a /b for which results are given is 0.8j this value is expected 
to lie close to the m aximu m that will be encountered on actual airplanes 
with fuselages . 


The setup of table I can of course be altered. For example, a 
solution may sometimes be required for the torsional stiffness necessary 
to maintain a given percentage of the rate of roll that would be obtained 
with a perfectly rigid wing. In such a case, the solution can be 
obtained. by recasting the last several columns, as shown in the table. 

If experimental data are available for any of the derivatives, these data 
may be used in place of the theoretical values of the present paper. 

Although an explicit solution for the reversal speed is impossible 
to obtain because of the complexity of the equations, the reversal speed 
is easily found by calculating the rate of roll through the speed range. 
However, the variation of the rate of roll with speed is usually of more 
interest than the reversal speed alone, so no extra calculations are 
involved in calculating the reversal speed. 


The twist parameter 0 r /& is directly proportional to q and 
inversely proportional to me r . (See equations (26).) Because of this 
relationship and because the other derivatives appearing in equations (2) 

and (3) are independent of q and mg , the values of 

, pb /2 V r 5 

and '^ p^2v) R found frQm table I for specified values of altitude and 

stiffness are easily converted to values for other altitudes and degrees 
of stiffness . 


Although the analysis was made by considering two wing panels lying 
in the same plane, experimental evidence exists which indicates that the 
results may with some justification be applied without change to configu- 
rations consisting of three or four wing panels, each containing an 
aileron, if the span and aspect ratio are taken as those of the wing 
formed by the reflection of a single panel. If the configuration consists 
of four panels, two of which contain ailerons, then the value of 00^ 

must be doubled. 
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The results of some calculations made for a hypothetical configu- 
ration are given in figure l 6 . 


CONCLUDING- EEMABKS 


An analysis has been presented of the effect of torsional flexibility 
on the rolling characteristics at supersonic speeds of tapered unswept 
wings with partial-span constant- percent- chord ailerons extending inboard 
from the wing tip. 

Because of the complexity of the equations resulting from the 
analysis, numerical calculations were made from the equations and the 
results aare presented in a series of figures. A computational form was 
prepared for use with the figures so that calculations could be made 
without reference to the analysis* 

Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

~ Langley Air Force Base, Va«, March 8 , 1949 
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Figure 7'.— Variation of rolling-moment— effectiveness coefficient with PA. — = 0.1 
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(From equation (20).) (k "blocks = 1 inch) 
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Figure 13.— Continued. (4 "blocks 
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Figure 13 •— Continued. (4 blocks = 1 inch) 




Figure 13*— Concluded. (4 blocks 
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Figure 15.- Concluded. (4 blocks = 1 inch) 






